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ABSTRACT 

We investigate the mid-infrared properties of the largest (42 objects) sample of radio-loud 
narrow-line Seyfert 1 (RL NLSl) collected to date, using data from the Wide-held Infrared 
Survey Explorer (WISE). We analyse the mid-IR colours of these objects and compare them 
to what is expected from different combinations of AGN and galaxy templates. We hnd that, 
in general, the host galaxy emission gives an important contribution to the observed mid- 
IR Hux in particular at the longest wavelengths (W3, at 12p.m, and W4, at 22pm). In about 
half of the sources (22 objects) we observe a very red mid-IR colour (W4-W3>2.5) that can 
be explained only using a starburst galaxy template (M82). Using the 22pm luminosities, 
corrected for the AGN conPibution, we have then estimated the star-formation rate for 20 
of these “red” RL NLSl, hnding values ranging from 10 to 500 Mq y '. Eor the RL NLSl 
showing bluer colours, instead, we cannot exclude the presence of a star-forming host galaxy 
although, on average, we expect a lower star-formation rate. Studying the radio (1.4 GHz) to 
mid-IR (22pm) flux ratios of the RL NLSl in the sample we found that in ~10 objects the 
star-forming activity could represent the most important component also at radio frequencies, 
in addition (or in alternative) to the relativistic jet. We conclude that both the mid-IR and the 
radio emission of RL NLS1 are a mixture of different components, including the relativistic 
jet, the dusty torus and an intense star-forming activity. 
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1 INTRODUCTION 

The importance of the nuclear active accretion phase (AGN) in the 
context of galaxy evolution has been full y recognized only in the 
last years (e.g see lHeckman & Bes3l2014l for a recent review). Ac¬ 
cording to the galaxy-AGN co-evolution scenario, at high redshifts, 
both star-forming (SF) activity and nuclear accretion must have 
proceeded simultaneously and at high rates thanks to the abun¬ 
dance of cold gas and to the frequent galaxy mergers that made 
the gas available for both processes. In a second phase, when the 
AGN luminosity was high enough to expel the gas from the galaxy, 
both star-formation and accretion rates gradually reduced even¬ 
tually leading to a “passive” elliptical hosting a quiescent super- 
massive black-hole. 

At low redshift, however, galaxy interactions and/or inter¬ 
nal galaxy dynamics can destabilize the residual gas and make 
it available for new epi s odes o f SF and nuclear activity at high 
accretion rates. iMathuil jlOOOh has proposed that the so-called 


Narrow-Line Seyfert 1 galaxies (NLSl), a peculiar sub-class of 
AGN that is characterized by rapidly accreting central black-holes 
with relatively small masses (between 10^ and 10* Mq ), may 
represent the low-redshift, low-luminosity analogues of the first 
high-z quasars, hosted by rejuvenated, gas-rich galaxies. In sup¬ 
port to this hypothesis, it has been noted that NLSl are more 
frequently associated to host galaxies with large levels of star- 
forming activity than Broad Line Seyfert 1 (BLSl) galax ies (e.g. 
IPeo, Crenshaw & KraemeJl2006l : ISani et aflbOloL l2012h . In this 
context, NLSl could be systems that are building up mass, evolving 
eventually into BLSL The time-scale of this process should depend 
on the radiative efficiency o f the accretion process, as discussed in 
lOrban de Xivrv et ^ ( 1201 ih . 


Interestingly, NLSl may be considered as “young systems” 
also from the radio point-of-view. NLS Is are typically radio¬ 
quiet (RQ) sources with only a small fraction that can be classi¬ 
fied as radio-loud (RL, ~7%, Komossa et aflllOObl and about 50 
in total known to date, iFoschini 2011 ). The radio properties of 
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these RL NLSl seem to be different from those observed in RL 
BLSl mainly because of the l ack of extended radio structures, 
except for very few cases (e.g. iDoi et al]|2012l) . It has been pro¬ 
posed that at least some of the RL NLSl can be associated to 
young radiogalaxies whose relativistic jet is still digging its way 
through the ISM of the host galaxy possibly evolving into a “classi¬ 
cal” giant radiogalaxy l Moranj 2000|^ Oshlack, Webster & Whiting 


200 ll: iGallo et aU l2006l: iKomossa et alj l2006l: I Yuan et alj 1200^: 


Caccianiga et alJl2014ll . Notably, young and compact radiogalax 


ies are observed also in the 3 radio-loud QSO discovered to date 
at z~6 jFrev et alj|201 ll . [2012ll . This further increases the possible 
analogies between local NLSl and high-z QSO. 

The possibility that NLSl are “young” systems where the 
AGN-galaxy co-evolution is taking place and where, in some cases, 
a young relativistic jet is ejected and observed in the early phase of 
its evolution, is intriguing: the study of these sources may offer a 
unique opportunity of investigating the co-evolution mechanism in 
both radio-loud and radio-quiet objects at relatively low-redshifts 
where the observational constraints are less challenging with re¬ 
spect to high-z QSO. 

In this paper we focus on the radio-loud fraction of NLS1 that 
is still poorly studied. In particular, we want to assess if RL NLSl 
are similar to RQ NLS 1 for what concerns their preference in be¬ 
ing hosted by star-forming galaxies. To date, only few RL NLSl 
have been studied from this point of view. Optical imaging has re¬ 
vealed the existence of arms or circumnuclear rings (possibly the 
consequence of a recent merger) and the presence of a circumnu- 
clear starburst in at least one RL NLS 1 (IH 0323-t342, Zhou et a^ 
1 20071 lAnton. Browne & March^l2008l: iLeon-Tavares et al.ll20liL 


In a recent paper ( Caccianiga et alj|2014h . a new case of RL NLS 1 

has been presented whose Spectral Energy Distribution (SED) is 
strongly suggestive of the presence of a quite intense star-formation 
(~50 Mq y“*). The presence of this intense SE was revealed by the 
mid-infrared (mid-IR) data at 12 and 22 q m from the Wide-hel d In¬ 
frared Survey Explorer (WISE) catalogue dWright et al.ll201(]h that 
showed a clear excess with respect to a typical dust emission from 
a radio-quiet AGN or to the IR non-thermal spectrum coming from 
a jet. We now want to extend such a result on a large sample of RL 
NLS 1 to establish whether the presence of star-formation is a pe¬ 
culiar characteristic of some isolated sources or, instead, is a global 
property of the entire class. 

The sample and the data used in this study are presented in 
Section 2. In Section 3 we discuss how the mid-IR colours (partic¬ 
ularly those including the 22 qm band) of the sources cannot be sat¬ 
isfactorily explained by the AGN emission alone. In Section 4 we 
consider the effect of galaxy emission on the observed IR colours 
and we infer the presence of a star-forming host in many objects 
of the sample. In Section 5 we compare the 1.4 GHz to 22qm flux 
ratios of the RL NLS 1 with those of blazars and infrared galaxies 
while in Section 6 we estimate the star-formation rate (SFR) present 
in the sources with the reddest W3-W4 colours. Conclusions are fi¬ 
nally presented in Section 7. 

Throughout the paper we assume a flat ACDM cosmology 
with Ho=71 km s“' Mpc^', Q\=0.7 and Qm^O.S. Spectral indices 


2 THE SAMPLE OF RL NLSl AND WISE DATA 

iFoschini et alj ( 1201 Sh has presented an extended study of the 



Figure 1. The WISE colours of the RL NLSl in the lFoschini et al.l l2015h 
sample. Open and filled circles indicate the RL NLSl with a redshift 
above and below 0.3 respectively. Sources circled by an open square are 
flagged as variable in the AllWISE catalogue while sources circled by tri¬ 
angles are flagged as extended. Different characteristic regions used to se¬ 
lect AGNs are also plotted: the blazar “WISE Gamma-ray Strip” (WGS) for 
BL Lacs and FSRQ (BZB and BZQ respe ctively, solid lines, gre en and blue 
in the electronic version), as defined by iMassaro et al.l i2012l) . the AGN 
wedge (dot-dashed line, magenta in the ele ctronic version) as defined by 
iMateos et al.l i2012l) and iMateos et al.l j2013h for X -ray selected AGNs, the 
mid-IR criterion (Wl-W2>0.8, brown dotted line) proposed bv IStern et alJ 
Finally, the red dashed line indicates the region covered by a simple 
power-law spectrum (slopes of Q'=0, 1 and 2 are indicated by tidemarks) 


42 RL NLSl with a large radio-loudness paramete£] (>10, see 
IFoschini et alJl20f^ for details) and a flat radio-spectrum (i.e. with 
blazar-like characteristics) discovered so fad Several pieces of in¬ 
formation, from the radio to the X-ray and gamma-ray band, have 
been collected and studied. The physical properties of the inner su- 
permassive black-hole (SMBH), like its mass and the accretion rate, 
have been also discussed. This is the largest sample of RL NLS 1 
collected so far and, therefore, it represents the best starting point to 
study the incidence of star-formation in the galaxies hosting these 
AGNs. It should be noted, however, that the type of selection of this 
sample, namely the requirement of a flat radio spectrum, strongly 
favours the presence of a relativistically beamed jet emission (as in 
blazars), which could be important also in the mid-IR band. This 


* The radio-loudness par ameter used in iFoschini et al.l l2015h is defined as 
the 5 GHz to 4400A flux density ratio. 

^ We note that only 22 sources out of 42 have a measured radio spectral 
index. Of these 22 objects, 20 have a flat slope (ctr 40.5) while 2 have 
a/i marginally above this value (0.55 and 0.58, respectively) but con sis- 
tent with it taking into acco unt the errors (see IFoschini et al.l 12015l and 
also iRichards & Listerll20fa for the recently measured spectral index of 
J0953+2836). The remaining 20 sources have no information on the radio 
spectral index and should be considered more properly as “flat-spectrum 
candidates”. 
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Table 1. WISE data on the sample of RL NLS1 


name 

other name 

z 

mag 

Wl 

(3.4|j.m) 

unc. 

S/N 

mag 

W2 

(4.6|j.m) 

unc. 

S/N 

mag 

W3 

(12 i+m) 

unc. 

S/N 

mag 

W4 

C22vim) 

unc. 

S/N 

JOlOO-0200 

FBQSJOlOO-0200 

0.227 

12.859 

0.023 

46.5 

11.754 

0.023 

46.7 

8.407 

0.024 

44.7 

6.187 

0.060 

18.0 

J0134-4258 

PMNJ0134-4258 

0.237 

12.157 

0.024 

46.1 

11.129 

0.021 

52.4 

8.125 

0.022 

49.1 

5.924 

0.042 

25.9 

J0324+3410 

1H0323+342 03 ‘ ^ 

0.061 

10.743 

0.022 

49.1 

9.791 

0.020 

54.7 

7.179 

0.016 

67.4 

4.847 

0.027 

40.1 

J0706+3901 

FBQSJ0706+3901 

0.086 

12.888 

0.024 

45.6 

11.959 

0.022 

48.9 

8.599 

0.026 

41.0 

6.049 

0.051 

21.4 

J0713+3820 

FBQSJ0713+3820 

0.123 

10.040 

0.022 

49.3 

8.990 

0.020 

53.5 

6.261 

0.015 

74.8 

3.984 

0.021 

51.7 

J0744+5149 

NVSSJ074402+514917 

0.460 

13.397 

0.024 

44.4 

12.384 

0.025 

43.7 

9.636 

0.044 

24.7 

6.829 

0.069 

15.6 

J0804+3853 

SDSSJ080409.23+385348.8 

0.211 

11.745 

0.023 

46.4 

10.750 

0.020 

54.3 

8.077 

0.022 

48.8 

5.372 

0.032 

33,5 

J08I4+5609 

SDSSJ081432.11+560956.6 

0.509 

14.221 

0.027 

40.4 

13.139 

0.028 

38.9 

10.411 

0.074 

14.6 

n.d. 



J0849+5108 

SDSSJ084957.97+510829.0 ‘ ^ 

0.584 

12.887 

0.024 

44.7 

11.956 

0.022 

48.8 

10.042 

0.049 

22.0 

7.552 

0.135 

8.0 

J0902+0443 

SDSSJ090227.16+044309.5 

0.532 

13.927 

0.026 

41.0 

13.096 

0.029 

37.7 

10.015 

0.064 

16.9 

7.058 

0.086 

12.7 

J0937+36I5 

SDSSJ093703.02+361537.1 

0.179 

12.165 

0.024 

45.7 

11.202 

0.021 

51.5 

8.135 

0.021 

50.8 

5.472 

0.035 

31.1 

J0945+1915 

SDSSJ094529.23+191548.8 

0.284 

12.042 

0.024 

45.9 

11.078 

0.022 

50.1 

8.279 

0.023 

47.4 

5.705 

0.044 

25.0 

J0948+0022 

SDSSJ094857.31+002225.4 - 

0.585 

13.282 

0.024 

44.5 

12.204 

0.023 

46.5 

9.096 

0.032 

33.4 

6.682 

0.088 

12.3 

J0953+2836 

SDSSJ095317.09+283601.5 

0.658 

14.809 

0.032 

33.9 

13.924 

0.039 

27.6 

11.335 

0.179 

6.1 

n.d. 



J1031+4234 

SDSSJ103123.73+423439.3 

0.376 

14.323 

0.029 

37.8 

13.300 

0.030 

35.8 

10.446 

0.078 

14.0 

8.392 

0.343 

3.2 

J1037+0036 

SDSSJ103727.45+003635.6 

0.595 

15.108 

0.037 

29.5 

13.854 

0.039 

27.8 

10.698 

0.103 

10.6 

n.d. 



J1038+4227 

SDSSJ103859.58+422742.2 

0.220 

12.549 

0.024 

46.0 

11.532 

0.021 

50.7 

8.870 

0.028 

39.1 

6.576 

0.058 

18.8 

J1047+4725 

SDSSJ104732.68+472532.0 

0.798 

14.749 

0.031 

35.2 

13.641 

0.033 

33.3 

10.435 

0.089 

12.2 

7.748 

0.189 

5.7 

J1048+2222 

SDSSJ104816.58+222239.0 

0.330 

13.499 

0.025 

44.3 

12.356 

0.026 

42.4 

9.292 

0.044 

24.5 

7.087 

0.114 

9.5 

J1102+2239 

SDSSJl 10223.38+223920.7 

0.453 

13.167 

0.024 

45.7 

12.042 

0.024 

44.7 

9.146 

0.034 

31.5 

6.489 

0.069 

15.7 

ill 10+3653 

SDSSJl 11005.03+365336,3 

0.630 

16.038 

0.058 

18.6 

15.230 

0.090 

12.1 

n.d. 



n.d. 



J1138+3653 

SDSSJl 13824.54+365327.1 

0.356 

14.055 

0.026 

42.1 

13.139 

0.028 

39.1 

10.532 

0.084 

13.0 

8.007 

0.207 

5.2 

J1146+3236 

SDSSJl 14654.28+323652.3 

0.465 

14.099 

0.027 

39.8 

13.208 

0.029 

38.0 

10.509 

0.080 

13.6 

8.571 

0.313 

3.5 

J1159+2838 

SDSSJ115917.32+283814.5 

0.210 

13.451 

0.025 

42.6 

12.343 

0.023 

46.3 

9.015 

0.034 

32.2 

5.982 

0.042 

25.8 

J1227+3214 

SDSSJ122749.14+321458.9 

0.137 

11.441 

0.023 

48.2 

10.328 

0.020 

53.7 

7.420 

0.018 

60.8 

4.802 

0.024 

44.8 

J1238+3942 

SDSSJ123852.12+394227.8 

0.623 

15.327 

0.037 

29.7 

14.525 

0.050 

21.9 

11.700 

0.183 

5.9 

n.d. 



J1246+0238 

SDSSJ124634.65+023809.0 

0.363 

14.163 

0.029 

37.4 

13.094 

0.029 

37.2 

10.346 

0.084 

12.9 

7,517 

0.161 

6.7 

J1333+4141 

SDSSJ133345.47+414127.7 

0.225 

13.102 

0.023 

47.2 

11.947 

0.022 

49.5 

8.763 

0.024 

44.6 

5.846 

0.043 

25.1 

J1346+3121 

SDSSJ134634.97+312133.7 

0.246 

13.744 

0.025 

43.1 

12.850 

0.026 

41.8 

9.693 

0.039 

27.6 

7.283 

0.105 

10.3 

J1348+2622 

SDSSJ134834.28+262205.9 

0.918 

14.710 

0.029 

37.0 

13.368 

0.028 

38.2 

10.400 

0.064 

16.9 

8.574 

0.287 

3.8 

J1358+2658 

SDSSJ135845.38+265808.5 

0.331 

13.133 

0.024 

44.7 

12.069 

0.022 

48.3 

9.225 

0.030 

36.2 

6.531 

0.058 

18.6 

J1421+2824 

SDSSJ142114.05+282452.8 

0.538 

13.674 

0.025 

43.0 

12.607 

0.024 

45.8 

9.614 

0.038 

28.7 

7.194 

0.096 

11.3 

J1505+0326 

SDSSJ150506.47+032630.8 

0.409 

14.028 

0.027 

40.1 

13.093 

0.029 

37.8 

9.840 

0.045 

24.4 

7.061 

0.086 

12.7 

J1548+3511 

SDSSJ154817.92+351128.0 

0.479 

13.770 

0.026 

42.1 

12.771 

0.024 

44.8 

9.883 

0.035 

31.4 

7.200 

0.072 

15.1 

J1612+4219 

SDSSJ161259.83+421940.3 

0.234 

13.308 

0.024 

45.8 

12.199 

0.022 

48.9 

8.135 

0.018 

62.0 

5,572 

0.030 

36.3 

J1629+4007 

SDSSJ162901.30+400759.9 

0.272 

13.290 

0.024 

45.6 

12.177 

0.022 

49.1 

9.471 

0.037 

29.4 

6.865 

0.070 

15.6 

J1633+4718 

SDSSJ163323.58+471858.9 ■ 

0.116 

12.135 

0.023 

46.6 

11.318 

0.020 

53.2 

7.642 

0.017 

64.0 

4.974 

0.024 

44.8 

J1634+4809 

SDSSJ163401.94+480940.2 

0.495 

14.883 

0.028 

39.0 

13.884 

0.031 

35.1 

10.711 

0.049 

22.1 

8.081 

0.119 

9.2 

J1644+2619 

SDSSJ164442.53+261913.2 

0.145 

13.283 

0.024 

44.7 

12.294 

0.024 

44.4 

9.455 

0.040 

27.1 

7.019 

0.100 

10.8 

J1709+2348 

SDSSJ170907.80+234837.7 

0.254 

13.758 

0.026 

42.5 

12.785 

0.024 

44.4 

9.903 

0.044 

24.6 

7.343 

0.117 

9.2 

J2007-4434 

PKS 2004-447 

0.240 

13.445 

0.025 

43.7 

12.349 

0.024 

44.8 

9.467 

0.035 

31.1 

7.025 

0.098 

11.0 

J2021-2235 

IRAS 20181-2244 ‘ 

0.185 

11.880 

0.023 

48.0 

10.997 

0.021 

51.6 

7.601 

0.017 

63.5 

4.710 

0.034 

32.2 


Notes: ^=flagged as extended in one or more bands in the AllWISE catalogue (ext_flg=3 or 5); ^=flagged as possible variable in W1 and W2 band in the AllWISE catalogue (var_flg>7) 


may partly (or totally) hide the emission from the host galaxy. 
To date, however, only few RL NLSl without “hla zar-like” char¬ 
acteris tics have been discovered and, therefore, the Iposchini et alj 
bOLSli sample is the only possibility for this type of analysis, at the 
moment. 

Th e mid-IR data have been collected from the WISE all-sky 
catalog jWright et alj|20l3) and, in particular, we have used the 
last version of the catalogue, the AllWISE data release (November 
2013), which has an enhanced photometric sensitivity and accuracy 
compared to the 2012 WISE All-Sky Data Release. All the NLS 1 of 
the sample are detected (S/N>3) in the WISE survey at 3.4 and 4.6 
micron (W1 and W2 bands, respectively), all but one are detected 
at 12p.m (W3) and 37 are detected at 22p.m (W4). Data for the 42 
RL NLSl are reported in Table [T] 

Four sources are flagged in the catalogue as possibly extended 
i.e. the source profile fit has a large (> 3) value in one or 
more bands. In one of these sources (J0324-I-3410/1H 0323-b342) 
imaging in R and B has revealed the existence of extended struc- 
tures (arms/rings) that could be th e consequence of a recent merger 
jAnton, Browne & March^l2008h . In the remaining 3 sources the 
available optical images (from SDSS or UKST) show the presence 
of a nearby (<5") object, possibly in interaction with the NLSL 
Finally, 3 sources in the catalogue (J0324-I-3410, J0849-b5108 and 


J0948-I-0022) are flagged as variable in the W1 and W2 band and 
for two of them (10849+5108 and J0948+0022) intra-day variabil¬ 
ity has been studied bv ijiang et al.l (l2012h . 


In Fig. [T] we show the distribution of the WISE colours in the 
Wl, W2 and W3 filters for the 41 RL NLSl detected in all the 3 
bands. In this figure we also show some characteristic regions re¬ 
cently adopted to select different types of AGN (see figure caption 
for details). Not surprisingly, all the objects fulfill the AGN defini¬ 
tion (Wl-W2>0.8) proposed bv IStem et alj ilOl'j) an d all but 2 lie 
in the AGN wedge discussed by Mateos et alj ( I20l3 . l2013h . Most 
(78%) of the RL NLS 1 occupy the region where blazars are usually 
found (blue and green solid regions) while 9 sources lie outside it. 
One source, in particular, shows a very blue W2-W3 colour (W2- 
W3=1.914), which is significantly different from the values ob¬ 
served in the other objects. This source, however, is one of the 3 ob¬ 
jects flagged as variable (J0849 +5108, one of the sou rces detected 
also in gamma-rays by FERMI. [Foschini^t^ l2015h. Indeed, the 
Wl and W2 magnitudes reported in Jiang et au i 2012h . based on 
the previous version of the WISE All-Sky catalogue, are signifi¬ 
cantly fainter (A ~ 1.1 mag) than the ones reported here. Therefore, 
the position on the W1-W2 vs W2-W3 plot of this objects is likely 
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Figure 2. The expected positions on the W1-W2 vs W2-W3 (left) and W1-W2 vs W3-W4 (right) plot of AGNs with different shapes of IR spectrum: 1) a 
simple power-law (black tick line) with slope ranging from 0.5 to 1.5; 2) a smooth broken power law, with slope before the break ranging from 0.5 to 1 and a 
slope of 1.5 after the break and with the break position varying within the WISE bands (dark blue shaded area); 3 ) a RQ AGN emissio n described by the 3 QSO 
templates from the SWIRE library (blue lines labelled as QSO, BQSO and TQSO, see text for a description and IPolletta et al.l2007l for more details) and with 
a redshift ranging from 0 (indicated by a cross) to 0.9 (the maximum value of the sample); 4) a combination of RQ AGN emission plus a power-law spectrum 
(light blue shaded area). In particular, we show the result of all the different combinations of slopes (from 0.5 to 1.5), redshift (from 0 to 0.9) and relative RQ 
AGN to power-law normalization. The arrows in the upper-right corners show the expected effect of an extinction of Av=2 mag at the mean redshift of the 
sample. An Av=2 mag is the maximum value expected for type 1 AGNs. We do not plot the an'ow along the W3-W4 since the extinction has a negligible 
effect in this case. Points as in Fig.[T] 


affected by the observed flux variations and, for this reason, it is not 
indicative of the real spectral distrihution of the sourc^ 


3 AGN EMISSION AND WISE COLOURS 

Even if the RL NLSl are broadly clustered around the region ex¬ 
pected for a power-law (PL) spectrum (with a spectral index Uir 
ranging from ~0.5 to ~1.5), many departures from this region are 
also observed: some sources occupy the region above this line (to¬ 
wards redder W1-W2 colours) while others occupy a region bluer 
in W1-W2 and redder in W2-W3. 

If the WISE emission of these RL NLSl were totally dom¬ 
inated by the non-thermal jet we would expect a power-law or a 
broken power-law spectral shape. In Eig. [2] (left) we show the re¬ 
gion on the W1-W2 vs W2-W3 plot where we expect to hnd a 
source whose IR spectrum is described by a broken power-law with 
a slope before the break ranging from 0.5 to 1 and a slope after the 
break of 1.5 and with the position of the break varying within the 
WISE bands (dark blue shaded area). As clear from Eig.[2(left), we 
expect that a source with a broken power-law spectrum would fall 
above the PL line. Several sources are actually found in this region 


although the majority of them seem to be much “redder”, in terms 
of W1-W2 colour, than what is expected from a broken power-law. 

In principle a significant level of extinction could change the 
position of a source with a power-law spectrum moving it away 
from the expected region, towards redder values of W1-W2 and 
W2-W3. However, since these sources are, by dehnition, type 1 
AGNs, i.e. objects where both the AGN continuum and the broad 
emission lines are visible in the optical spectrum, the extinction is 
not exp ected to have an im portant impact on the mid-IR colours 
(e.g. see lMateos et alj|2015h . In the upper-right side of both panels 
of Pig. 1^ we indicate with an arrow the expected etfect of an extinc¬ 
tion of Av=2 mag at the mean redshift of the sample (<z>=0.36). 
An Ay of 2 magnitudes is th e maximum expected for type 1 AGNs 
as discussed, for instance, in ICaccianiga et all ( 1200 Sh . 

It is also possible that part (or the totality) of the ob¬ 
served IR emission is produced by the radio-quiet compo¬ 
nent of the source (i.e. the dusty torus) and not by the non- 
thermal jet. Indeed, the presence of this component is ex¬ 
pected and it has been revealed in m an y radio-loud AGN 
and blazars (e.g. Malmrose et al.l 20n Raiteri_et^ 201A|; 


iGurkan, Hardcastle & Jarvis! |2014| ; ICastignani & De Zotti l2015h . 
In Pi^ 2](left ) we plot the position of 3 radio-quiet type 1 QSO tem- 
plateU (from IPolletta et al]|2007h in the range of redshift between 


^ It should be considered that the AllWISE catalogue, that we are using 
in this work, has been produced by combining all the data from both the 
WISE cryogenic and post-cryogenic survey phases, the latter being carried 
out only in the W1 and W2 bands. Therefore, the periods over which the 
images have been collected and combined are different for the 4 bands 


'* These templates have been derived bv IPolletta et al.l feoOTl) by combin¬ 
ing the SDSS quasar composite spectrum and rest-frame IR data of a sample 
of optically-selected type 1 QSOs observed in the SWIRE program. The 3 
templates have the same optical spectrum but three different IR SEDs: a 
mean IR spectrum, obtained from the average fluxes of all measurements 
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0 and 0.9 (the range of z of the sample, red solid lines). The re¬ 
gions occupied by these 3 templates lie above the PL line and reach 
redder W1-W2 values than the broken power-law, in particular for 
high redshift sources (z between 0.5 and 0.9). Notably, the source 
with the reddest value of W1-W2 (=1.34, J1348+2622) is also the 
object with the largest redshift in the sample (z=0.918). Only the 
TQSOl template can reach such a red W1-W2 colour that a bro¬ 
ken power-law, instead, does not seem able to reproduce. We con¬ 
sider this as an indication that, at least in a number of objects, the 
mid-IR emission is probably mainly produced by the radio-quiet 
component (i.e. the torus) rather than the relativistic jet. 

Finally, we consider also the possibility that both the dusty- 
torus and the jet emissions are present in the spectrum. To this end, 
we combine the 3 QSO templates with a power-law emission, at 
different values of redshift (from 0 to 0.9), varying the relative nor¬ 
malizations of the two components (light blue shaded area). 

From the analysis of Fig. (left) we conclude that the sources 
lying above the PL line in the W1-W2 vs W2-W3 plot could be 
explained either by a broken power-law emission, expected when 
the jet emission is dominant, or by a RQ AGN (i.e. torus) emis¬ 
sion or by a combination of the two components. Only one object 
(J08494-5108) lies above the PL but in a isolated region toward 
blue W2-W3 values (Wl-W2=0.931 and W2-W3=L914). As said 
before, this is a highly variable source and, therefore, its position is 
not indicative of the actual IR spectral shape. 

In ~ 10 cases, the objects fall well below the PL line and, there¬ 
fore, their colours are not easily explained by a torus/jet emission 
nor by a combination of the two. The number of sources not ex¬ 
plained by these models is even higher if we consider the longest 
wavelength magnitude (W4, at A=22[im) in the diagnostic plot i.e. 
if we plot W1-W2 vs W3-W4 (Fig. [2] right). As before, we over¬ 
plot the expected regions followed by a power-law, a broken power- 
law and a RQ AGN component plus a power-law emission. From 
this figure it is clear that only few objects can be explained by these 
models while the large majority (~ 80-90%) lie well below the PL 
line thus requiring a different (very red) component. As described 
in the following section, the best candidate for the emission of this 
red component is the host galaxy. 


4 THE EFFECT OF GALAXY DILUTION ON IR 
COLOURS 

We consider here the effect of the combination of an AGN plus the 
host galaxy emission on the WISE colour plots. To th is end, we use 
some representative templates of the SWIRE library jPolletta et alj 
l2007h Specifically, we use an early-type galaxy (“E115” template), 
an Sc galaxy (“Sc” template) and a starburst galaxy (“M82” tem¬ 
plate, Eig. [3- We start from a galaxy template and add a PL with 
different slopes and varying the relative AGN/galaxy normaliza¬ 
tion. In Fig.|4]we show the predicted effect of the host galaxy di¬ 
lution on a power-law emission. An early-type host is expected to 
move the AGN colours towards the left/bottom side of the plot i.e. 
it is expected to produce bluer colours in both W1-W2 and W2-W3 
(W3-W4) values as it gets more and more important. A starburst 
host, instead, is expected to make the W1-W2 values bluer than 


(“QSOl”), a template with high IR/optical flux ratio template, obtained 
from the highest 25 per cent measurements per bin (“TQSOl”), and a low- 
IR emission SED obtained from the lowest 25 per cent measurements per 
bin (“BQSOl”). 


W1W2 W3W4 



X(//m) 

Figure 3. Templates from the SWIRE library jPolletta et al.ll2007l) used in 
the analysis. On the top axis we indicate the centres of the four WISE bands 
at the average redshift of the sample (0.36). 


the power-law while both the W2-W3 and the W3-W4 colours be¬ 
come redder, thus bringing the object on the right-most side of the 
plot. Other starburst templates show a similar behaviour. This is 
the effect of the strong IR emission of the starburst component that 
becomes more and more important in the W3 (and W4) bands. 

In Fig. |4] we also show the position of 
SDSSJ143244.91 - b30143 5.3, a RL NLSl discussed in 
ICaccianiga et alj ( l2014h whose SED has revealed the pres¬ 
ence of a star-forming host galaxy with a SFR of 50 M© y“'. In 
order to correctly reproduce the infrared SED (in particular the 
W4 photometric poi nt) of this object we used the M82 template by 
IPolletta et alj JtOO^ . 

Fig. m (left) shows that the host galaxy may be contributing 
significantly to the mid-IR emission in a number (10) of sources 
of the sample, i.e. those objects that occupy the region below the 
PL line. For the sources just below this line the contribution from 
the host galaxy in these WISE bands is not enough to distinguish 
the type of the host galaxy. In ~4 cases, instead, the very red W2- 
W3 colours seem to favour a star-forming host galaxy. We note 
that the position of SDSSJ143244.91-I-301435.3 on this plot does 
n ot clearly indicate a star-f orming host galaxy. Indeed, as discussed 
in ICaccianiga et all ( l20I4h . the strongest indication for a M82-like 
host galaxy comes from the W4 band (22 qm) which shows a clear 
excess. Fig. |4] (right) shows that using the W4 filter we can better 
discriminate between the different types of host galaxy. In partic¬ 
ular, the mid-IR colours of the 22 sources with W3-W4>2.5 (like 
SDSSJ143244.91-H301435.3) can be explained only by combining 
the AGN emission to a starburst galaxy host. 

In order to better show the effect of the host galaxy on the 
observed spectrum we report in Fig. |5] the average UV-to-mid-IR 
SED of the NLS1 of the sample divided according to the value of 
W3-W4 colour (greater than 2.5 or lower than 2.3, in order to keep 
more separated the two groups). The average SEDs have been cre- 
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Figure 4. The effect of the host galaxy emission on a power-law spectrum, in the WISE colours plots (W1-W2 vs W2-W3, left, and vs W3-W4, right). The 3 
regions between the long-dashed/short-dashed/point -dashed lines repres ent the expected paths of a combination of power-law plus a galaxy template taken from 
the SWIRE library (M82, Sc and E115 respectively, IPolletta et alJl2QQ7h . varying the relative intensity of the two components, the power-law slopes (ranging 
from 0.5 to 1.5) and the redshift of the source (±i(T=±0.2 from the mean value of z). Note that the limits are different from those of Fig.[2|in order to better show 
the different paths produced by different host galaxies . Symbols and arrows as in Fig. [2 The blue star indicates the position of SDSSJ143244.91+301435.3, a 
RL NLSl studied in detail bv ICaccianiga et aU i2014h and whose SED is suggestive of the presence of a starburst host galaxy. 




Figure 5. Average UV to mid-IR SED of the RL NLSl in the sample divided according to the W3—W4 colour. Red points represent the average values of 
Log/lL,i in each bin of LogT while the shaded area represents the region including 68 % of points. The thick black line is a model composed by a reddened 
(the Ay value is indicated in the label) AGN template (QSOI from IPolletta et al.l2Q07h plus a galaxy template: in the left panel an early-type galaxy template 
(E115) is used while, in the right panel, we adopt a starburst galaxy template (M82). 
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ated by normalizing each rest-frame SEDs at 4ttm and computing 
the average value of Log/lL,i in bins of Logd. The grey area re¬ 
ports the region including 68% of points. The da ta used to produce 
these SED are taken from iFoschini et alj ( 1201 Sh and include, be¬ 
sides WISE photometry, data from Swift/UVOT observations an d 
from exi sting public surveys like SDSS fDR9. lAhn et alJ|2012h . 
2MASS dCutri et alj|2003h and USNO-B dMonet et alfcOOak 

To correctly interpret the shape of these average SEDs we 
model them with the sum of a galaxy and a QSO template, vary¬ 
ing the relative normaliza tion, the amount o f nuclear extinction and 
the host morphology (see lBallo et alJ[20l^ for a more detailed de¬ 
scription of the procedure). We consider different combinations of 
extinctions and galaxy/AGN templates and we find that templates 
with a low IR emission (early-type/Sc galaxies) are able to describe 
the observed average SED of the objects with bluer W3-W4 colours 
(left panel) while these templates systematically fail in reproducing 
the W3/W4 data points in the average SED of sources with redder 
W3-W4 (right panel). In this latter case, instead, a starburst galaxy 
template, like M82, is more suited to reproduce the steep increase 
of the data points towards longer wavelengths. 

We note that the two average SEDs show some differences 
also in the UV/optical part. We can explain this difference with a 
different (average) extinction level. While to model the average- 
SED of the sources with a blue W3-W4 colour we have to apply 
only a small level of extinction (Av/=0.1 mag), in the case of the 
red W3-W4 objects, we need a higher value of Ay (0.35 mag). The 
explanation of this result could be related to the enhanced presence 
of dust in highly star-for ming galaxies wit h re spect to less star¬ 
formi ng galaxies (see e.g. iBuat et al.|[2003 and iDommguez et alj 
l2013h . 


5 JET OR STAR-FORMATION? 

So far, we have focused our analysis on the UV to mid-IR data 
and, in particular, on the mid-IR colours. From these data we have 
inferred the presence of a star-forming/starburst host galaxy in at 
least half of the RL NLSl in the sample. Since an intense SF is ex¬ 
pected to produce an important emission also in the radio band it is 
possible that the radio luminosities observed in these sources are in 
part contaminated by the SF activity. In principle, it is even possible 
that a high SFR can explain the totality of the observed radio emis¬ 
sion. This would be an unexpected result, since the sample under 
investigation is composed by objects with large radio-loudness pa¬ 
rameters and, in some cases, a flat radio spectrum. These properties 
are usually considered as strong indications that the radio emission 
is likely produced by a relativistic jet pointing towards the observer 
(and, thus, relativistically beamed). 

Since jet-dominated objects and star-forming galaxies have 
significantly different radio-to-mid-IR flux ratios, it is interesting 
to compare this parameter computed for the RL NLS1 to the values 
usually observed in jet-dominated AGN (i.e. blazars), from the one 
hand, and to the values observed in star-forming galaxies, on the 
other hand. 

To quantify the radio-to-mid-IR flux ratios we use two equiv¬ 
alent parameters. The two point spectral-index, defined between 
1.4 GHz and 22p.m (rest-fram^ as: 


^ For the K-coiTection we use, for each source, the slope computed between 
12 and 22|j.m (observed frame) fluxes and we assume Q'=0 in the radio band. 
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q22 

1 0 - 1-2 



Figure 6. Distribution of the two-point spectral index between radio 
(1.4 GHz, rest-frame) and mid-IR (22qm, rest-frame) of the RL NLSl (blue 
histogram) and of blazars (BZQ) in the BZCAT (black shaded histogram). 
This latter histogram has been re-scaled by a factor 10 to facilitate the com¬ 
parison. The flux densities have been K-corrected assuming a slope of 0 in 
the radio and using the observed IR slope between 12 and 22 qm. The filled 
blue histogram shows the RL NLS 1 detected in the gamma-ray by FERMI. 
The dashed red histogram includes only the red (W3-W4>2.5) objects. The 
top axis reports the values of q22. For comparison, the typical range of q24 
(defined at 24qm instead of 22qm) observed i n the infrared galaxies is also 
indicated (light blue region. lRieke et al.l2009h . 


^-og{S xACBzISiivim) 

®14 =- - ( 1 ) 

Eog(V|.40H;/''22qm) 

which is typically used in blazars studies, and the parameter 
q22 defined as: 

q22 = Log(S 22 qm/‘S i. 4 G«:) (2) 

The radio-to-mid-IR indices for the 37 RL NLSl in the sample 
detected in the W4 band are reported in Fig.[3 For comparison, we 
also report the distribution of o-j^ for a sample of 430 blazars with 
a quasar spectru m (Flat Spectrum Rad io Quasar, FSRQ) extracted 
from the BZCAT dMassaro et al.ll200^ and detected with high sig¬ 
nificance (S/N>10) in all the WISE filters. To help the comparison, 
this latter histogram has been re-normalized by a factor 10. We also 
indicate the typical range of r adio-to-mid- I R valu es observed in the 
infrared galaxies reported by iRieke et alj l l2009l) . In particular, we 
report the observed range of q24, a quantity similar to q22 but de¬ 
fined at 24qm instead of 22qm using Spitzer MIPS data. The aver¬ 
age value found is 1.2-1.4, depending on the infrared luminosity of 
the source, with a dispersion of 0.24 dex. 

The first consideration on Fig.is that the RL NLSl span a 
wide interval of radio-to-mid-IR indices, covering the entire range 
between the values observed in infrared galaxies to those observed 
in blazars. This is a first evidence that the radio and mid-IR emis¬ 
sion of these sources are a mixture of different components. In 
sources with steep indices (>0.2), i.e. within the range of values 
observed in blazars, the relativistic jet is likely the dominant com- 
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ponent in the radio band (while in the mid-IR different emission 
could be present, including the dusty torus and/or the host galaxy). 
In sources with <-0.25 (i.e. q22> 1) the SF activity may be con¬ 
tributing significantly to the observed radio emission. This is likely 
true, for instance, in the sources with red W3-W4 colours (dashed 
red histogram), that are the majority in this interval of q22 values. 
The most obvious interpretation for the objects showing intermedi¬ 
ate values, instead, is a combination of jet and SF (plus the dusty 
torus, in the mid-IR). 

In support to the hypothesis that the relativistic jet is domi¬ 
nant in sources with steep values we note that 6 out of the 15 
RL NLSl with af-^ >0 (i.e. 40%) are detected in gamma-rays by 
FERMI (filled histogram in Fig. [IJ, something that is usually con¬ 
sidered as an indication for the presence (and the dominance) of 
a relativistic jet, while only one object (J1102-1-2239) out 22 with 
4 <0 (5%) is a FERMI source. The dominance of the relativistic 
jet in 4 of the FERMI detected RL NLS 1 with q22>0 (J0324+3410, 
J0849-I-5108, J0948-I-0022 and J1505-I-0326) has been recently 
quantified also by lAngelakis et alj l l2015h by means of a system¬ 
atic monitoring in the radio and millimetric bands. 

This confirms the presence and the importance of the relativis¬ 
tic jet in the sources with a positive value of while it seems to 
suggest that other mechanisms are at work in addition (or instead 
of) to the relativistic jet in most of the remaining sources. It is worth 
noting, however, that a low value of (or a high value of q22) 
does not necessarily exclude the presence of a jet as demonstrated 
by the detection in gamma-rays by FERMI of one of the objects 
with a high value of q22 (=1.34, J11021- 22390 Another exam¬ 
ple is the sources J1227-1-3214 that, in spite of a l arge q22 value 
(=1.2 9), has an inverted radio spectrum (ap=-1.04. lEoschini et alj 
l2015h something that suggests a non-thermal origin of the radio 
emission (or at least of a fraction of it). It should be considered that 
variability can play an important role in determining the properties 
of a source at different wavelengths, including its position within 
the histogram of Fig.[^ during high activity periods, a source may 
appear as jet-dominated in many (if not all) frequencies, showing 
for instance a flat radio spectrum, a strong gamma-ray emission or a 
high radio-to-mid-IR flux ratio, while, during low activity periods, 
the signs for the presence of the jet can be hidden by the RQ AGN 
emission or even by the host galaxy light. For all these reasons, si¬ 
multaneous multiwavelengths monitoring campaigns, together with 
radio follow-ups, in particular at VLBI resolution, are fundamental 
to disentangle the different components and establish their relative 
importance at different wavelengths. 


6 ESTIMATE OF THE STAR-FORMATION RATE 

In this section we want to quantify the intensity of the star- 
formation at least in those sources where the mid-IR colours are 
suggestive of the presence of a starburst host. Specifically, we con¬ 
sider here only the 22 sources with a red W3-W4 colour (W3- 
W4>2.5). To limit the contamination due to the AGN nucleus (ei¬ 
ther the dusty torus or the jet component) we analyse here the 
longest wavelength WISE band (W4 at 22pm) since, according to 
the results previously discussed, this is the band where the star¬ 
forming nature of the host galaxy seems to emerge more clearly. 
We first compute the 22pm luminosities (AL^). We derive the 

® Interestingly, this source was indicated bv iFoschini et al.l l2015h as pos¬ 
sible outlier in the gamma-ray/radio correlation. 



Figure 7. Star-formation rate estimated from the 22pm luminosity cor¬ 
rected for the AGN contribution (see text for details). This estimate has 
been done for the objects whose WISE colours suggest the presence of a 
star-forming host galaxy (W3-W4>2.5). Shaded histogram includes only 
the sources with a q22> 1. 

22pm flux densities from the W4 magnitude^ and applied a K- 
correction using the spectral slope computed between 12 and 22 pm 
(observed frame) fluxes. The computed 22pm luminosities of the 
sources range between 10*** to lO'^ L© with a large fraction (76%) 
of objects with L 22 pm >10" Lq. 

In order to derive an estima te of the SFR in th ese sources we 
can use the relation presented in iRieke et alj ( l2009l) between SER 
and L 24 pm (their eq.lO and eq.ll). However, since the L 22 pm lu¬ 
minosity is likely contaminated by the AGN, in particular the dusty 
torus component, we have to apply a correction. To estimate the 
importance of the AGN emission at 22pm we proceed as follows. 

Jet contribution. We estimate the contribution from the jet 
starting from the radio flux density measured at 1.4 GHz, assuming 
the average value of the observed in the blazars (ESRQ) of the 
BZCAT catalogue {af;^=0.27, see Fig. [6) and subtract it from the 
observed flux density at 22pm (rest-frame). For the 2 objects that 
have >0.27 the estimate of 22pm due to the galaxy is not pos¬ 
sible and we exclude them from the analysis. For the sources with 
radio-to-mid-IR flux ratios similar to those observed in IR galax¬ 
ies (q-i^^ <-0.25), the correction for the jet emission is negligible 
(<1%) and has no impact on the final results. As discussed in the 
previous section, in these sources even the radio emission could be 
mainly (or totally) attributed to the star-formation activity. 

Torus contribution. We estimate the contribution from the 
torus using the flux densities at 4.6pm (rest fram0) and assuming 

^ For this conversion we assume an IR spectral slope of 1. The differences 
in the derived flux densities when assuming significantly different slopes, 
from -1 to 2, are below 0.7% 

* For the K-correction we use here, for each source, the observed slope 
computed between 3.4 and 4.6 pm. 
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name Log SFR 

(Moy-*) 

(1) (2) 

Log SFR unc. 
(M0y-‘) 
(3) 

q22 

(4) 

W1-W2 

(5) 

W3-W4 

(6) 

O'R 

(7) 

107064-3901 

0.95 

[ 0.92 . 

, 0.98 ] 

0.82 

0.93 

2.55 

- 

J0744-r5149 

2.50 

[ 2.43 . 

, 2.56 ] 

0.47 

1.01 

2.81 

- 

J0804+3853 

1.95 

[ 1.94. 

. 1.96] 

1.48 

0.99 

2.70 

- 

J09024-0443 

2.70 

[ 1.70. 

,3.70] 

-0.65 

0.83 

2.96 

0.07 ± 0.01 

J0937+3615 

1.92 

[ 1.91 . 

, 1.93 ] 

1.31 

0.96 

2.66 

- 

J09454-1915 

2.14 

[ 2.06. 

, 2.21 ] 

0.60 

0.96 

2.57 

- 

J1102+2239 

2.48 

[ 2.47 . 

, 2.49 ] 

1.34 

1.13 

2.66 

- 

J1138+3653 

1.16 

[0.16. 

,2.16] 

-0.15 

0.92 

2.52 

0.5 ± 0.09 

J1159+2838 

2.11 

[2.10. 

,2.11 ] 

1.39 

1.11 

3.03 

- 

J1227+3214 

1.73 

[ 1.71 . 

, 1.74] 

1.29 

1.11 

2.62 

-1.04 + 0.07 

J1246+0238 

1.83 

[0.83 . 

, 2.83 ] 

-0.36 

1.07 

2.83 

0.55 ± 0.06 

J1333+4141 

2.17 

[2.16. 

,2.17] 

1.39 

1.16 

2.92 

- 

J1358+2658 

2.11 

[ 2.09 . 

,2.12] 

1.30 

1.06 

2.69 

- 

J1548+3511 

2.18 

[ 1.18 . 

,3.18] 

-0.77 

1.00 

2.68 

0.26 ± 0.01 

J1612+4219 

2.35 

[ 2.35 . 

, 2.36 ] 

1.32 

1.11 

2.56 

- 

J1629+4007 

1.39 

[ 1.06. 

, 1.72] 

0.29 

1.11 

2.61 

-0.68 ± 0.02 

J1633+4718 

1.83 

[ 1.74. 

, 1.92] 

0.23 

0.82 

2.67 

0.42 ± 0.01 

J1634+4809 

2.07 

[ 1.93 . 

, 2.20 ] 

0.15 

1.00 

2.63 

- 

J1709+2348 

1.36 

[ 1.33 . 

, 1.39] 

0.96 

0.97 

2.56 

- 

J2021-2235 

2.48 

[ 2.46. 

, 2.50 ] 

0.81 

0.88 

2.89 

0.5 ± 0.07 


Table 2. Estimates of the star-formation rate for the sources with W3-W4>2.5 and with < 0.27, based on the 22|j.m luminosities corrected for the AGN 
contribution (see text for details). Col.l: source name; col.2: logarithm of the SFR derived from the 22|j.m luminosity; col.3: uncertainty interval on the SFR 
that takes into account only the errors on the estimate of the 22|j.m luminosities of the host galaxies (see text for a description). In addition to this uncertainty 
it should be considered also the one related to the SFR-F(24|j.m) relation (about ~0.2 dex. iRieke et al.l 2009|); col. 4: logarithm of the 22|j.m-to-1.4GHz flux 


density ratio; col. 5: W1-W2 colour; col. 6: W3-W4 colour; col. 7: radio spectral slope (Sy cc v “) from 
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that this is entirely produced by the AGN (jet and torus). This is a 
reasonable assumption since, as discussed in the previous sections, 
in the W1 and W2 filters we are probably preferentially observing 
the AGN emission (see also the average SEDs of Fi g. [5). This is 
also co nsistent with what has been recently found by Mateos et alj 
ll2015h studying the mid-IR properties of a sample of X-ray selected 
AGN. If also the host galaxy is contributing to the 4.6p.m emission 
then the final estimate of the 22p.m emission from the host galaxy, 
and hence of the SFR, would be even higher. As for the 22[im flux 
density, also the 4.6[J.m flux density has been corrected for the jet 
contribution, using the average 1.4 GHz-to-4.6p.m spectral index 
observed in the FSRQ of the BZCAT (< ttj ® >=0.43). Using the 
resulting flux density and the slope between 4.6 and 22p.m mea¬ 
sured on the QSOl template from SWIRE (q'= 1.1), we can then 
estimate the contribution of the torus at 22[im and subtract also 
this component from the observed flux density. 


We find that the average contribution from the torus at 22p.m 
is more important (~40%) than that from the jet (~7%). In order 
to take into account the uncertainties related to the initial assump¬ 
tion on the of the jet and the IR slope of the torus, we have 
repeated these computations with different input values: aj^=0.15 
and 0.39 (i.e. the Icr interval observed in the FSRQ of the BZCAT) 
and a;R=1.0-l.l, for the infrared slope of the torus, corresponding 
to the extreme values measured in the 3 QSOl template of SWIRE 
(TQSOl, BQSOl and QSOl). We take the extreme values obtained 
in all these combinations as the uncertainty interval of the 22p.m lu¬ 
minosity from the host galaxy. The computed 22p.m luminosities of 
the host galaxies range from 10*** to 5xlO"F0. 

We then use the 22p.m luminosity corrected for the AGN 
co ntribution to estima te the SFR using the relation presented 
in iRieke et ^ ( l2009l) . their eq.lO and eq.ll, and assuming 
L24(j.m(M/P5) ~L 22 p.m- The results are reported in Tab.|5]and plot¬ 


ted in Fig.|7] The SFR of these 20 objects are above ~10 Mq y“' 
and, in several cases, larger than 100 Mq y“', similarly to what is 
usually observed in Luminous Infrared Galax ies (LIRG) or even 
in Ul traluminous Infrared Galaxies (ULIRG, I Sanders & Mirab^ 
Il996h . This result confirms and quantifies what we have inferred 
from the analysis of the mid-IR colours, i.e. the presence of an in¬ 
tense star-forming activity in at least half of the RL NFS 1 in sam¬ 
ple. 

We finally note that a SFR estimate for the RL NLSl with 
W3-W4<2.5 would not be reliable since the evidence of a signif¬ 
icant contribution of a star-forming host galaxy to the emission at 
22 qm is not clear. Running the analysis described above on the 
blue objects, we find that, on average, only a small fraction (~8%) 
of the observed 22 qm luminosities can be due to the host galaxy 
the remaining fraction being produced by the torus (56%) and by 
the jet (36%). The estimate of the 22p.m luminosity due to the host 
galaxy (and therefore the SFR quantification) depends too critically 
on the starting assumptions on the of the jet and on the IR spec¬ 
tral shape of the dusty torus to give reliable results. The only stable 
quantities that we can derive are the upper limits on the SFR based 
on the total observed 22p.m luminosities. These values do not ex¬ 
clude the presence of an intense SFR (> 10 Mq y“*) also in the 
objects with W3-W4<2.5. From these upper limits alone it is not 
possible to understand whether the host galaxies of these AGN are 
intrinsically dilferent from those observed in “red” RL NLSl or, 
rather, their SF activity is simply hidden by a particularly strong 
nuclear light. We have, however, other pieces of information that 
can help to distinguish between these two hypotheses. 

In Fig. [5] we have presented the average SEDs of red (W3- 
W4>2.5) and blue (W3-W4<2.3) objects, that clearly show differ¬ 
ent shapes at A > 10p.m. While the average SED of the red objects 
requires the presence of an M82-like template to be correctly re- 
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produced, the average SED of the blue objects can be modelled 
using an Sc or E115 template. An M82 template would over-predict 
the emission at A >10pm. We note that the normalization of the 
host galaxy template is quite well constrained in the spectral region 
around ~-lpm (rest-frame). This suggests that the blue objects are 
hosted, on average, by different types of host galaxy. 

In order to assess whether there are differences also in the lu- 
minosites of the AGN (either the radio-quiet or the jet component) 
between red and blue objects we can compare the average luminosi¬ 
ties at 4p.m, where the torus emission is expected to dominate, and 
in the radio band. At 4p.m the average values are only marginally 
different, with the blue objects being about 30% more luminous 
than the red ones. Even in the radio band we do not have evidence 
that the jets in the blue objects are more luminous than the jets 
present in the red sources. Rather, the average 1.4 GHz power of 
blue objects is a factor ~1.6 fainter than in red objects. The aver¬ 
age 22 qm luminosity, instead, is significantly different bewteeen 
the two groups of sources, where the red objects are almost a factor 
2 more luminous that the blue sources. 

Overall, these numbers seem to indicate that the lack of a clear 
detection of the host galaxy in the blue objects is not due to a more 
luminous radio jet and/or AGN (radio-quiet) but, rather, to a less 
luminous 22qm emission, i.e. a less important star-forming compo¬ 
nent at work in these sources. As explained above, the upper limits 
on the SER are in any case large so we cannot exclude the pres¬ 
ence of an intense SE also in these sources but, on average, their 
activity, if present, seems to be less intense than the one observed 
in the host galaxies of the red RL NLS1. This picture needs a con¬ 
firmation through an accurate follow-up (in particular of the blue 
RL NLSl) in the mid-IR or in the far-IR spectral region, where the 
maximum output of the SE activity is expected. 


7 CONCLUSIONS 

We have analysed the mid-IR colours of a sample of 42 RL NLS 1 
using WISE data with the aim of quantifying, for the first time, the 
star-forming activity in the galaxies hosting this class of AGNs. The 
sources cover the redshift interval from 0.06 to 0.92. Thirty-seven 
objects are detected in all the 4 WISE bands. In order to under¬ 
stand the origin of the mid-IR emission in these sources we have 
compared their positions on different diagnostic plots with those 
expected from a combination of AGN and host galaxy templates. 
The results can be summarized as follows: 

• In general, the RL NLS 1 of the sample occupy a region in the 
W1-W2 vs W2-W3 plot (commonly used in the literature) that is 
typical of emission line AGN, both radio-quiet and radio-loud. 

• WISE colours can be reproduced by a combination of AGN 
and host galaxy emission, the relative importance of the two com¬ 
ponents depending strongly on the band under consideration. While 
the emission from 3.4qm up to 12qm (observed frame) of most 
(3/4) of the sources can be explained by the AGN alone (either the 
relativistic jet, the dusty torus or a combination of the two), at 22qm 
the emission usually requires the (red) contribution from the host 
galaxy. 

• Twenty-two sources show very red W3-W4 colours (>2.5) that 
can be explained only by assuming the presence of a young and 
“active” host galaxy, i.e. a starburst galaxy similar to M82. In the 
remaining cases, the intensity of the AGN emission in the mid- 
IR bands does not allow to unambiguously assess the type of host 
galaxy. 


• The 22qm luminosities of the NLSl detected in W4 are typ¬ 
ically above lO" Lq. We have used these luminosities, corrected 
for the AGN contribution, to estimate the SER for 20 out of the 
22 sources with W3-W4>2.5. The computed values of SER range 
between 10 and 500 Mq y“'. These values are similar to thos e ob¬ 
served in LIRG or even in ULIRG iSanders & Mirabeilll996h . Eor 
the sources with W3-W4<2.5 the presence of a star-forming host 
galaxy cannot be excluded although the SER is expected to be, on 
average, lower than in red sources. 

• Although the sample studied here has been designed to in¬ 
clude NLS 1 with a relativistic jet (possibly pointing towards Earth), 
by means of a selection of sources with a large radio-loudness pa¬ 
rameter (radio-loudness para meter >10) and, pref erentially, a flat 
radio spectrum (a^ <0.5, see lEoschini et alj|2015l for details), our 
analysis shows that in the mid-IR the jet is not necessarily the domi¬ 
nant component, even in sources that have been detected in gamma- 
rays by EERMI. In particular, we have estimated that, on average, 
at 22qm only ~20% of the luminosity of the sources with W3- 
W4>2.5 comes from the jet, the remaining part being emitted by 
the dusty torus and by the host galaxy. 

• Even at radio frequencies (1.4 GHz) the emission of a number 
(~10) of sources in the sample, i.e. those with high mid-IR-to-radio 
luminosity ratios, is not necessarily due to the jet but it is likely 
produced, in part or entirely, by the SE activity. 

In conclusion, our analysis has shown that RL NLS 1 are of¬ 
ten associated to “active” host galaxies, with SER in the typical 
range of LIRG/ULIRG sources. Therefore, from this point of view, 
RL NLSl are more similar to RQ NLSl than to radio-loud AGNs 
with broader emission lines, usually hosted by “passive” elliptical 
galaxies. These results support the idea that NLSl, both RQ and 
RL, are systems in the early phase of their evolution, when the host 
galaxy is experiencing a high level of star-forming activity and, at 
the same time, the central super-massive BH is rapidly accreting 
and building up mass. Studying the properties of these sources and 
their differences with respect to the other radio-loud AGNs will 
help us to understand the role of relativistic jets in the AGN/galaxy 
evolution and to unveil any possible form of feedback (either posi¬ 
tive or negative). To this end, simultaneous follow-ups, both in the 
mid-IR and in the radio, are mandatory in order to accurately disen¬ 
tangle the different components at work in these active nuclei and 
to study their possible interplay. 
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